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\l ' Abstract. The radial equation of a simple potential model has long been known to yield 

an exponential survival law in lowest order (Breit-Wigner) approximation. We demon- 
strate that if the calculation is extended to fourth order the survival law exhibits a 
parabolic short time behavior which leads to the quantum Zeno effect. This model has 
further been studied numerically to characterize the extra exponential time parameter 
which compliments the lifetime. We also investigate the inverse Zeno effect. 

dl! PACS. 03.65.Xp 



^ : I. INTRODUCTION 

The exponential survival law is known to be an excellent phenomenological fit to unstable phenomena. 
However, there is no rigorous derivation of this law in quantum mechanics. Most text book derivations 
KJi ' are "classical" in nature since they refer only to decay probabilities. In quantum mechanics, one would 

, require an exponential time dependence for the survival amplitude [1] . To achieve this one must make 

' approximations. Standard procedure is to consider a tunnelling process, as we do in this paper. Some 

earlier important theoretical papers upon tunnelling are given in Refs. [2]. 

Quantum mechanics allows us to say that the survival probability P{t) is definitely not exponential 
for short and long times. For short times it can be argued that a power expansion in t must lack 
the linear term, i.e. {dP/dt)t^„ = 0. This result when combined with the hypothesis of "frequent" 
measurements leads to what is known as the quantum Zeno effect (QZE) [3]. It was first named 
the quantum Zeno paradox by Misra and Sudarshan [4] precisely because it was considered a false 
result [5]. Nowadays, the QZE is generally acknowledged as a real phenomenon and indeed there are 
a number of experiments which claim to have verified it and others are planned [6]. Once accepted, it 
is even possible to predict a QZE within a classical calculation [7]. As for the long time behavior, this 
is predicted to be a power law behavior in t . This latter result can be derived even for a Breit-Wigner 
spectrum by imposing a low energy cut-off [8] (which is expected on physical grounds). The long-time 
property of P(t) will not be treated specifically in this paper but it is referred to in the concluding 
Sections. 

Returning to the short time behavior, the simplest demonstration of the non-exponential behavior 
of the survival law is based upon the hermitian nature of the Hamiltonian [9] (when all channels are 
considered). Consider a state \tp{t)) initially in a (quasi-bound) state represented by |'0(O)) = \ipo)- It 
is assumed that \ipo) can "decay" into another or other states. At time t the state will have evolved 
into 

1^(0) = exp{-iHt) IV'o) , 
where H is the Hamiltonian. Expanding in a power series of t, 



\m) = [I - tHt - ^ + 0{f)] iVo) 
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The amplitude for non-decay is given, up to a possible phase, by 

where ( )o stands for the average over the state |V'o)- Whence, the non-decay probability is, 

= l-z{{H- H+))o - i i{H')„ + {H+')„ - 2 {H)„{H^)„) f + 0{f) . 
Now, using H = wc find 

The linear term in t which corresponds to the linear term in H has cancelled. Of course this derivation 
assumes the existence of {H^)o , in addition to (-ff)o- It is useful to observe at this point that if this 
demonstration where extended to all powers in t then not only would the linear term vanish but all 
odd powers of t would vanish. Wc shall return to this when we discuss some numerical calculations 
in Section IV. 

The QZE merits a name (even if this choice is not quite appropriate) because it implies a potentially 
spectacular phenomenon - the prediction that " frequent" short time tests of the state of an unstable 
system will inhibit its decay. Much has been written upon this, in particular upon the question of 
what may constitute an "observation" or test of the system. We wish to add here only a comment 
upon this fascinating subject. A measurement of a system normally produces a collapse of the wave 
function. It is therefore perfectly reasonable to expect that observations should modify, for example, 
the survival law. However, the exponential curve e"*/"^ is unique in this respect because it is not 
altered by measurements. This fact is connected to the mathematical feature that the average value 
{{t — t„)) from t„ to oo (lifetime) is independent of the lower limit. This is the reason one does not 
need to know when the, otherwise identical, unstable particles in a "sample" where created in order 
to measure their lifetime. Even if each had been created at a different time, one can treat them as if 
"newly created" at the conventional time tg. In more colorful terms, a series of operations of "cut and 
paste" (measurements) are undetectable only if performed upon a single exponential function. The 
QZE is thus a consequence of a particular example of the non exponential nature of P{t). What is 
somewhat peculiar is that the absence of the linear term in t is more important (because in principle 
more easily verifiable) than the absence of any other or indeed of all the other odd powers of t. 

In this paper we wish to examine the survival law with the aid of a particular potential model 
(Section II), already used in the literature for a derivation of the exponential law, albeit as an ap- 
proximation [10]. In Section III, we will repeat this derivation with particular emphasis upon the 
assumptions/approximations that lead to the exponential law. We will then go beyond the lowest or- 
der calculation (Breit-Wigner form) to a two pole approximation. Recently a two pole approximation 
in a quantum field model was introduced by Facchi and Pascazio [11]. In Section IV, we shall show that 
at this more sophisticated level we predict the necessary conditions for the QZE. Indeed, we speculate 
in the conclusions that only the lowest order approximation loses this effect. In Section V, we present 
the results of some numerical calculations and in particular the dependence of an extra time constant 
upon the parameters of the model. We will phenomenologically parameterize these dependencies. We 
conclude in Section VI with the resume of our results, a discussion of the "inverse Zeno effect" [12] 
(described in Section V) and some observations upon the significance of the extra time parameters. 



II. THE POTENTIAL MODEL 

The starting point of our analysis is the three dimensional Schrodinger equation for a particle of mass 
TO in a spherically symmetric potential V{r), i.e. a function only of the magnitude r of r. Explicitly, 
we will use 
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V{r) 



Region 1 



Region 2 



Region 3 



V{t) = 



< r < a , 
a < r < b , 
b < r . 



a b 

While our subject matter is not a stationary state problem, it can be conveniently analyzed with the 
use of the stationary state solutions, each of which has a simple (phase) time dependence. Due to the 
spherical symmetry of the potential energy and its time-independence, the Schrodinger equation 



2m 



V{r) 



i;{r,t) 



can be separated by using the energy eigenstates 



9 ^ 



(1) 



where Y, 



equation 



arc the spherical harmonics. Thus Eq.(Q yields the well known time-independent 
1 / 5" 2 d\ + 



2m \ dr^ 



dr 



2 m 



Vir) 



(2) 



For a given value of E and angular momentum, there are two linearly independent solutions of the 
above second order equation, which at the origin go like or 1/r' + \ However, those which behave 
like l/r' + ^ must be rejected since it can be shown that ^ ^ is not a solution of the above 

eigenvalue equation for r = 0. This is because the Laplacian of YJ"*(0, 0)/r' + ^ involves the l-th 
derivative of the Dirac delta function 6{r) [1]. 

Henceforth, we shall limit our attention to s-wave solutions (/ — 0) and consequently, the previous 
equation reduces to 



1 

2m 



2 d_ 

r dr 



Vir) 



RE{r)^ER^{r) 



(3) 



with Re{Q) which, in accordance with our previous discussion, is a constant. If we now put 

Rsir) = UE{r)/r , 

we obtain the equation for the modified radial wave function Ue (r) 

u'^{r) =2m[Vir)~E]uEir) , (4) 

where the prime stands for d/dr. The acceptable solution of Eq.Q must go to zero at the origin, 
Me(0) = 0. The standard procedure for finding the stationary solutions can now be applied. The 
general solution for u^ir) and any E < Vq in the three regions are 



UEir) 



Region 1 
Region 2 
Region 3 



< 7' < a , 
a < r < b , 
b <r , 



Ci sin(A:r) 



[ k = \/2mE ] , 



ikr 



-ikr 



[p^ ^2miV„-E)] 



(5) 



The requirement of continuity of Ue and u'^ at r = a, 6 (a direct consequence of the differential 
equation for Ue) can be conveniently expressed in matrix form as 



Ci sin(fca) 
k Ci cos(fca) 



peP"" 



e-P" \ 




pe-P" J 


D2 



and 



^ikb 

ike''''' 



-ike-''''' 



D, 



pe 



pb 



e-p''\ 


'a' 


pe-P'') 
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Eliminating C'2 and D2, we obtain C3 and in terms of Ci, 



-,—ikh 
'ikh 



-ikh 



pe 



,pb 
pb 



-pe 



-pb 



,pa 



1 pa 



sin(fca) 
k cos(fca) 



(6) 



^From this equation it is straightforward to show that jCaj = I-D3I. The solutions in region 3 involve 
waves "moving" in both directions. 

The normalization of the ^^(r) is determined by the standard condition for a continuous spectrum 



dr w^(rX,(r) = 6{E - E') 



(7) 



We recall that the definition of Re{t) = Usir^jr eliminates the factor r'^ in the integration variables 
(d^r = r^drdf2) while the 47r from the angular integration is accounted for in the definition of the 
spherical harmonic 1^°. The above normalization condition is dominated by region 3 since it has 
infinite range and determines uniquely |C3|(= ID3I) 



(ICal' + lZJan / dr cy.vVi}i~'i^)r\^ mb{k~k') Ik 



{k,k' > 0) 



whence 



\D, 



m/ (2TTk) 



(8) 



With this condition we can calculate all the other coefficient moduli. We shall need in what follows 
the explicit value of |Ci|^. From Eq.®, after some manipulations, we find 



where, 



IC3I = i I cxp[p (6 -a)]a^{l- i f ) + exp[-p (b - a)]a_ {l + i ^) \ \C,\ 
a± — sin(fca) ± - cos(fca) . 



(9) 



Whence, by using Eqs. JHJ and © 

\C,\' = ^ { [cxp{2 apw)al + cxp{- 2 apw) at] (fc' + p') + 2a+a_ (k' ^ p') 



(10) 



where 

w = [b — a)/a . 

Now, exp(2ajOw) can be considered to be a very large number for the purposes of this study. This 
can always be guaranteed for E < Vq hy increasing the width 6 — a of the potential. Under these 
conditions, in general the energy eigenstates have small values of |Ci| because of the large value of 
the first term in the square brackets above. There exist however, quasi-stationary states defined by 
a+ = when the opposite is true. These quasi-stationary energy levels coincide with the discrete 
spectrum of an associated potential model, 
V{r) 



Region 1 


Region 2 







V(r) = 







< r < a , 
a < r . 



The energy eigenvalues of this auxiliary potential, indicated generically by E^ in the following, are 
given by the solutions of the transcendental equation = 0, i.e. 



tan(aA:o) = - fco/po 



(11) 



where fco = V 2 mEa and po = ^j2m{Vo — E,,). The explicit values for a given Va and well size a can be 
calculated numerically when needed. We shall name the eigenfunctions of this second potential Uo(?') 



( \ — J Region 1: < r < a , Ci sin(fcor) . 
""^"■^ ~ 1 Region 2: 



a < r , 



4 exp{-por) . 



(12) 
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The continuity conditions for Uo[r) and u[^{r) yield both the transcendental equation and 



Ci = exp(-poa) a/ kl + pi / ko 
Finally, from the normalization of Uo{r) we obtain 

|ci|' = 2po/{l + apo) . 



(13) 
(14) 



III. THE BREIT-WIGNER APPROXIMATION 

In the previous section, wc outlined the potential model and introduced the variables and wave 
functions, including that of the auxiliary potential with bound states Uo{r) whose energy eigenvalues 
equal those of the quasi-stationary energies of our potential. We begin here by assuming that initially 
our radial state is Uo{r). Since this is not an cigenstatc of our potential but of the auxiliary potential, 
it will not remain localized permanently but eventually leak into region 3 through tunneling. To be 
more precise Uo{r) already has a small tail outside the potential barrier, however this can be safely 
neglected as we do in Eo. fH^ below. 

We now decompose Uo{r, t = 0) into energy eigenstates UE{r) 



Mo(r, i = 0) = Uo(r) 



d^; goiE) u^ir) 



from which it follows that 



Uo{r, t) 



dE g„{E)uE{r)e- 



-lEt 



(15) 



The spectrum go{E) is explicitly derived in Eqs. H18|l and (|19() below . At any time t > 0, the amplitude 
for finding the particle still in the quasi-stationary state Uo{r) is given by 



(16) 



dr ul{r,0)uoir,t) , 
and hence the non-decay or survival probability P{t) is 

^dE g„{E)e-'''' I dE' gl{E') I dr <,(r)w^(r) 

Jo Jo 

which, after using the normahzation condition Eq.Q, becomes 



Pit) 



dE \g,{E)\ 



2 -iEt 



(17) 



The index upon ga{E) reminds us that we have chosen one of the possible many quasi-stationary 
states Eg. To determine go{E) we use 



dr<(r)uo(r,0)= / dE' g,{E') / dr ul{r)uE.{r) = g,{E) . 



The left hand side of this equation can easily be calculated if one assumes 6 » a so that Mo(r, 0) can 
be considered negligible for r > b. This is the first of the approximations made. Now, as we shall 
see, for suitable choices of parameters, go{E) is highly peaked around the quasi-energy Eo (second 
approximation compatible with the first). In these cases, the function u^ir) does not, practically, 
differ from Ua (r, 0) in region 1 and 2 except for its normalization. We may then approximate the 
above integral to an integral over only region 1 and 2 and use 



u*(r) w CiU*(r,0)/ci 



{r<b), 
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whence for E around En 



g,{E)^{CJc,) dr ul{r,Q)u,{r,0)^{CJc,) / dr <(r, 0) Mo(r, 0) = Ci/ci 



(18) 



otherwise 



9.{E) ~ . (19) 
By using Eqs. (|10|l and (|14|) . we then obtain (for E around Eg) 

\go{E)\^ ^ """i^t"""^ {[exp(2ap«;)<+exp(-2apH«-] (F + p^) + 2 a+a_ (P - p^) } . 

For convenience, wc now make a change of variables from E and to a and t , by defining the 
dimensionless quantities 

a = ak and t = t/(2ma^) . 

Henceforth, wc will indicate by Eo the chosen initial quasi-stationary eigenvalue. In these new vari- 
ables, 



where 



Pit) - [4(l + apo)/a7rpo]' 



da I cxp(— icT^t ) //(cr) 



(20) 



/(^) 



exp 



2 ^y2maFVo~a^ luj 



a cos (T 



exp ( —2 ■\/2 ma^VJ, — w ) ( sin ct — 



a~ COS" (T 
2 ma^Vo — (7^ 



V2ma^K - cr^ 

(T COS (T 



sm" (T — 



V2ma2yo - cr^ 



2maW„ 



(21) 



In accordance with our second approximation, we will assume that l/f{a) is peaked about (Tq corre- 
sponding to the quasi-stationary energy Eo {a„ = a\J2 niEo). We then expand /(cr) around do and 
keep only up to quadratic terms in Act = ct — CTq. The validity of the peaked form is not in question. 
It can eventually be verified a posteriori. However, our assumption neglects terms of order (Act)'^ and 
higher. We shall indicate this fact by a superscript 2 on /(ct), 



/f^l (ct) = exp(2 apo w)[{l + ap,) {K + pi) f [{o - o.f - e (ct - ct,,) + 7= ] /a^K p^ 



(22) 



where 



e = exp(— 2 apa w) 



{l + apo){k,-+p,')' 



and 



1 ^0 + Po 

' 2 1,2 _ „2 

'''0 Po 



By convention we choose 7 to be positive. The linear term in Act with coefficient e can be eliminated 
by a shift in variables yielding a Breit-Wigner form for the a (and hence energy) spectrum. However, 
we shall not do this here since we intend to generalize the above formula to fourth order in the next 
Section, and then it is not in general possible to eliminate by a shift both of the odd powers in Act. 

The quadratic expression in square brackets in Ea. (|22|l has two complex conjugate zeros which 
appear as poles in our integrand for P{t). These are at 



CT = CTo 



Po 



ko + Po 

^0 i ^ Uo 1 



7 ± i 



2 fcoPo 

K + pS 



(23) 



with Xq and Do the real and imaginary parts of ct respectively. The integral in Eq. H2()|l is thus propor- 
tional to 

n + oo 

dCTexp(— iCT^t) [ (ct — Zo)(ct — z*) , 
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with Zo = Xo+iyo- Since by assumption this integral is sharply peaked about CTo we can formally extend 
the lower limit of the integral from to — oo. This is the third approximation in this derivation. The 
theorem of residues then yields 

dz exp(~izH) [{z - Zo){z - z*)]"' = {t^/Vo) exp{-i z*^i) . 

Taking the square of the modulus of this we finally obtain, after some simple algebraic manipulations, 

PW^expf-^t) . (24) 



The above expression is of the desired exponential form 

P{t)=cxp{~t/To) , (25) 
normalized (automatically) to P(0) = 1, and with predicted lifetime 

(26) 



2 ^0^0 

The narrow width of the Breit-Wigner assumed in the derivation is guaranteed if 7 <C 1 which 
corresponds to exp(2 apo w) 3> 1. This is just the condition anticipated in the previous Section. Thus, 
for the validity of the above derivation the choices of i?o, Vo, and 6 — a must be such as to satisfy this 
condition. In all our subsequent numerical calculations this is indeed the case. 



IV. HIGHER ORDER CORRECTIONS 

The calculation in the previous Section is a non-relativistic quantum mechanical derivation of the 
exponential law. It is incompatible with the existence of an experimental QZE and if repeated for 
negative (as well as positive) values of t it yields the non-analytic expression, 

P(t)=cxp[-|t|/r] . 

This satisfies the symmetry property P{—t) = P{t) (time reversal invariancc) consistent with the 
absence of all odd powers of t (Section I), but it cannot be expressed as a power series in t. 

What is surprising, at first, is the following: If P{t) is calculated by evaluating the integral in 
Ea. (|20|l numerically for various i, and then, for small t, fit with a polynomial, we find the odd powers 
of t greatly suppressed. That is, the numerical calculation is far nearer to what one may have expected, 
a priori, on the basis of the arguments in our introduction. This P„u„,(i) is practically an even function 
in t, and in particular the linear term in t is almost zero. To understand this we recall the expression 
of the non-decay amplitude, see Ea. H17(l . 



£dE \g„{E)\-' e-^^* « £dE \go{E)\' [cos(M) + z sin(M)] 



the terms in square bracket are bounded and by assumptions the |(7o(-E')|^ is highly peaked around 
E = Eg and tends to zero at least like {E — Eq)^^ when E ±00. The odd terms in t lie in the 
sine term which is pure imaginary. Hence, on taking the square of the modulus, all odd powers of t 
disappear and our numerical calculation simply reflects this. Indeed the presence of the very small 
odd terms in t can simply be attributed to rounding errors. Note that this result is derived for a Breit- 
Wigner spectrum, which is consequently not in itself a sufficient condition for obtaining the exponential 
law. Thus our numerical integral of the Breit-Wigner would seem to be a better approximation to the 
curve P(t) than the analytically derived exponential law. The following higher order approximation 
confirms how difficult it is to derive a simple exponential result. 

We now generalize the calculation of the previous Section by expanding /(u) up to fourth order 
in Act . We evaluate [a) as an improved approximation to the Breit-Wigner case. The rest of the 
procedure remains the same. Up to second order we obtained Ea. (|24|l . Now, we must evaluate 

dz exp{—i z'^t) [{z — Zi){z — z*){z — Z2){z — z*)] ^ , 
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where in our subsequent plots the values of Zi^2, with real and imaginary parts Xi^2 and j/i, 2 respectively, 
will be obtained numerically. However, even without these calculations we expect, from the second 
order approximation, that one of the poles must correspond to Zg. Let this be Zi {z^ ~ Zg). The other 
then provides a deviation from the exponential law (small if 2/2 3> j/i). Indeed analytically 



dz 



exp(— I z^i ) 



Defining 
we find 



{z ~ Zi){z ~ z*){z - Z2){z ~ z*) z* 
cxp{ia) = (z* - Z2)/izi - z*) 



exp{-i z*'i) exp(-iz*^i) 
4 I Viizt - Z2) y2{zi-zl) 

and 13 ^xl-xl + yl-yl , 



P^^Ht) = N I exp(-4xiyi<) + — exp(-4 2:22/20 + 2 — cxp[-2 (xij/i + 2:22/2) t] cos(a + /3?) 
where the normalization at t = 0, P(0) = 1, fixes 



(27) 



vl 



2 — cos a 

V2 



From Ea. l|27(l we see that as yi/2/2 the single exponential is indeed recovered. This justifies 
our prediction that Zi =2:01 in agreement also with our numerical derivations of z^. Apart from the 
coefficients in front of the additional terms in P''''(i) we observe that for 2:22/2 ^ x-^y^ (valid in all our 
numerical solutions for z-^ and Z2) the exponentials in these terms fall off far more rapidly than the 
leading exponential in t. Thus, they are of interest only for small t. 

The surprise here (at least for the authors) is that P'*' {t) displays exactly the absence of a linear 
term, notwithstanding that it is still an approximation, only a step above the Brcit-Wigncr approxi- 
mation. This is independent of the specific expressions for Zi_2 which would be far too complicated to 
be given algebraically in terms of the parameters. To see this, we expand P'*'(t) as a power series in t. 
The coefficient of the linear term t after a little algebra can be shown to be zero identically. To avoid 
any misunderstandings, wc emphasize that this fact has nothing to do with numerical integrations. 
We have used the theorem of residues to evaluate the integrals in the above. Anticipating one of our 
conclusions, one therefore sees that one must make a very precise set of approximations to derive the 
exponential law. 



V. NUMERICAL POLE CALCULATIONS AND INVERSE ZENO EFFECT 

So far we have not derived any numerical values for the poles. In this Section, we shall do so with the 
principle objective of finding a phenomenological fit to one of the time constants that characterizes 
the improvement to the exponential curve. In the expression for P'*'(t) three parameters appear x-^y-i, 
2:22/2, and (3 . The lifetime t is no longer sufficient by itself to parameterize P{t). One would like to 
interpret these constants, or more precisely, the related time parameters, in physical terms. 

In general the pole positions are complicated functions of the input parameters a, b, VJ, and Eg of 
the model. The only one we give here is an approximate algebraic expression for the P'^'(t) lifetime 
(to) derived in Section III. By using Eas. l|23|l and H26|) . wc find 



kaPo 
^0 ~ Po 



( 2 (To e + ) 



Dropping the term proportional to e 

2 



16 [EgiVg^Eg) 



,3/2 



l + a^2m{Vo- Eg) 



exp 



2{b~a)^2m{Vg~Eo) 



(28) 



(29) 



To study the other time parameters it is far easier to select the input values and solve numerically 
the equation for the pole positions. 

P'\<J)=0. (30) 
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It is convenient, as far as possible, to work with dimensionless quantities. In Table 1, we show for 
various values of Vo the solution of the above equation for and z^. The multiple of tt (a convenient 
choice suggested by the transcendental equation) that appears in the header of each sub-table corre- 
sponds to the number of quasi-stationary energies. For each value of these energies we list, for selected 
w, the values of Xi^a, !/i,2, and the dimensionless parameters Ti^^/ma^ defined below. The values of z^ 
coincide within numerical accuracy with those of Zo, the pole position for /'^'(cr) = 0. Consequently, 
to a good approximation, this pole still determines the lifetime for the process 



T, = ^ . (31) 



The other pole produces the modification from the simple exponential curve and determines a second 
exponential time parameter <^ {x^yi <C X2y-2) 

= 7y ■ (32) 

In Fig. 1 we display the dependence of Ta upon w for various Eg and an arbitrary chosen value of Vq. 
In Fig. 2 we show the dependence of Ta upon ^J2miy^^^^I^ for different values of w. From the above 
curves we see that, to a very good approximation, grows linearly with w and hence with the barrier 
width b— a. It is also proportional approximately to [ 2 m{Vo — Ea) \~^''^ . This for the parameter ranges 
considered. We defer an interpretation of these results to our conclusions. 

For t ^ 0, the exponential survival law must necessarily lie below the physical curve P{t) if this 
lacks a linear term since both are normalized to P(0) = 1. The fact that the lifetime measures the 
mean time for decay means that an exponential curve with said lifetime as input must eventually rise 
above P{t) to compensate for the small time "short- fall", i.e. there is a cross-over point. As an aside, 
we note that the power law in t predicted for P{t) at long times means that any exponentially falling 
curve must lie below P(t) asymptotically adding another cross-over point. The first cross-over allows 
for the so-called inverse Zeno effect (IZE) [12]. This says that if a first time measurement is made 
after this cross-over point (but logically before any second cross-over point) the system will have been 
found to decay at a rate greater than " expected" . 

First, let us analyze critically this definition. Unlike for the QZE this measurement or these mea- 
surements (if repeated at the appropriate time intervals) after the inverse Zeno cross-over only contra- 
dicts our expectations based upon the use of an exponential curve with the experimental lifetime. But 
the choice of exponential is somewhat arbitrary. We can give at least two alternative proposals for the 
reference exponential curve. The first, which is the most natural phenomenologically, is to refer to an 
exponential best fit to P{t). This fit almost certainly will not have the exact same lifetime as P{t). 
The second possibility is that one may compare P{t) with a theoretical single exponential (possibly 
from a model calculation as in our case). These possibilities distinguish themselves from the original 
because they do not necessarily have a cross-over point and hence need not imply the conditions for 
an IZE. 

We have a natural choice of exponential in our model, the Breit-Wigner exponential. We have thus 
looked for cross-over points by confronting P''''(t) with P'^'(i). In all cases examined they have been 
found. Thus, even with our choice of exponential the IZE is possible. We can see this more clearly by 
observing that, in our tables, yi/2/2 ^ 1 and Xi ~ x^, thus we can approximate Eq. ()27|l as follows 



P''"(t) w exp(-4xiyit) | l-f- 2 — [eyiY>{-2x2yJ^) cos{a + (it) - cosa]^ 
Consequently, 



(33) 



where 



« 1 + 2 ^ ^{t) , (34) 

I N / x^y^ \ ( P \ 

ip{t) ~ exp —2 t cos a H t — cos a . 

^ ma'' / \ ma' I 



Intersection points occur when Lp{t) = 0, i.e. when 

/3 



exp (2 t] = cos { a^ t] I cos a . (35) 

V ma' J \ ma' I 
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The trivial solution is, of course, t = 0. The condition for the IZE can be determined from the relative 
values of the time derivative of the functions on the left and right hand sides of the above equation 
at t ~ 0. The necessary and sufficient condition for a non-trivial intersection is 



2x22/2 < — P tana . 



(36) 



This is because, with this condition satisfied, the exponential starts off with a flatter growth than the 
oscillating and bounded right hand side of Ea. (|35|) . and hence they must necessarily cross. 

If Xi = X2 sin a = 0), we find only the trivial solution t = Q. On the other hand, from an 
examination of Table 1 , we see that 



■ (3 tan a = 2 



+ yl-vl){xi- X2){yi+V2) 



{xi - X2) - ivi + ^2) 



2 2^22/2 



1 - 






\ ■2' 1 -3^2 / 



> 2 X2y2 



Hence a cross-over occurs for all our tabulated cases. In Fig. 3 we show examples of this by plotting 
ip{t) against t for several quasi-stationary states of a fixed barrier height Vq and an arbitrary chosen 
value of w. The large dots identify the cross-over times. 



VI. CONCLUSIONS 

We have discussed in this paper the survival law and its relevance to the QZE and the IZE within the 
context of a simple potential model. We have shown that the general arguments that allow for the 
QZE based upon the hermiticity of the Hamiltonian are indeed valid in this model unless one makes 
very specific assumptions and/or approximations. In particular, the exponential law, which does not 
allow for the QZE, is obtained only if two approximations are made: 

1 - The denominator /(cr) in the spectral integral for the amplitude of non-decay is so highly peaked 
about the initial quasi-state that only the lowest - second order - terms in Act need be considered; 

2 - The integral in energy E or equivalently a is extended below the threshold to —00. 

Point 2 is made in order to apply the theorem of residues essential for deriving a single potential 
survival law. It is also the cause of loss of analyticity in t, although this is normally ignored by 
arguing that only positive values of t are physically significant. 

Our analysis has gone on to show that while point 1 (the Breit-Wigner spectrum) is not the sole 
cause for the theoretical lose of the QZE, it is the primary one. Indeed, if one includes higher order 
corrections, such as the fourth order terms in Act, then the resultant P''''(t) displays the absence of 
linear term in t, however small the extra contributions might be. 

This result was somewhat of a surprise. One could have reasonably supposed that only in the limit 
of an exact calculation, i.e. when all orders of Act are included, is the linear term absent. Indeed the 
authors only expected that to fourth order the coefficient of the linear term would be reduced when 
compared to that of the exponential curve. Instead we have been able to show that it is rigorously null 
already in P'*i(t). Note that to prove this we have applied points 2 and 3. Point 2 is also significant 
for another reason. It is known to be the cause for the lose of the long-time power law behavior, 
and this is independent of any other approximation made. Our conclusions upon the question of the 
theoretical prerequisite for the QZE is that it is indeed a feature of P{t) and that one must make very 
specific assumptions or approximations to avoid this effect. Only with the explicit exclusion of short 
and long times, can a single exponential curve be a good approximation to P{t). 

We have also considered in the previous Section the question of the existence or otherwise of the 
IZE. There is in our opinion a measure of ambiguity in the condition for the existence of this effect 
because of a non unique choice of the exponential used as reference. Since any exponential is at best 
an approximation to the physical curve, its definition is subject to discussion. An exponential with 
the same lifetime as the real curve is the conventional choice. However, there are other choices for 
which the possible existence of an IZE is far from obvious. We have shown that one possibility, the 
comparison of P'*''(t) with P'^'(t), docs indeed allow for an IZE. Nevertheless, we do not consider the 
IZE and the QZE comparable phenomena. The former (IZE) refers to our "expectations", and as we 
have argued our expectations are subject to some ambiguity. The latter (QZE) predicts a physical 
suppression of decay through continuous observation and it is completely independent of the existence 
or otherwise of an approximate exponential curve. 
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Of particular interest in our results is the dependence of the second exponential time parameter . 
The first time parameter is essentially the lifetime of the system and has a characteristic exponential 
rise with barrier width 6— a, whereas this second time parameter can be phenomcnological represented 

by 

ma^ w m (b — a) 

T2 ~ — , = — , . (37) 

^2ma-{Vo-Eo) ^2m{Vo-Eo) 

This form is very suggestive. Consider an energy eigenstate above the barrier, £' > V^. In the barrier 
region {a < r < b) the particle will have velocity v = ^2ni{E — Vo)/ni. Consequently the time taken 
in crossing the barrier will be: 

mlb — a) , ^ 

At = , ^ ' =. 38 

This "mirrors" the expression for T2 (except for the feature that Eg has a discrete spectrum). Both 
tend to infinity if the particle energy and potential energy are set equal. A type of consistency con- 
dition. Furthermore, both grow linearly with barrier width. This tempts us to speculate (and at this 
stage it is no more than an hypothesis) that T2 is a measure of the time that the particle takes to 
tunnel free. Equivalently, we may say that T2 is the time during which the particle is neither bound 
nor free [7]. However, since this infringes upon the question of transit times and super-luminal veloc- 
ities [13], which is a very topical, but also complex subject, we desist from any further considerations 
at this point. The tau are not the only time parameters in P'*'(<), we also have 1//3 in the oscillatory 
term. Mathematically it is the natural consequence of interference between the two exponentials of 
the non-decay amplitude. It is therefore a wave-like property of the particle. 
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Table 1. Numerical values of xi,2 , J/1,2 and ri,2 for chosen values of Vo- Each block lists the parameters for 
all the quasi-stationary eigenvalues Eq. 
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Fig. 1. Dependence of T2 on w for several of the quasi-stationary energies Eo of a fixed Vo 
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Fig. 2. T2 vs l/\/Vo — Eg for various quasi-energies Eg and w values. The lines drawn demonstrate excellent 
linear fits, for a fixed w. 




Fig. 3. The function ip{t) vs t for the indicated values of Vo and a;. The five curves are for different values of 
Eq . For each of the curves the cross-over points are indicated by the large dots. 



